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Abstract 

We analyze asymptotically a differential-difference equation, that arises in a Markov- 
modulated ffuid model. Here there are identical sources that turn on and off, and 
when on they generate fluid at unit rate into a buffer, which process the fluid at a rate 
c < A^. In the steady state limit, the joint probability distribution of the buffer content 
and the number of active sources satisfies a system of -|- 1 ODEs, that can also be 
viewed as a differential-difference equation analogous to a backward/forward parabolic 
PDE. We use singular perturbation methods to analyze the problem for N ^ oo, with 
appropriate scalings of the two state variables. In particular, the ray method and 
asymptotic matching are used. 

1 Introduction 

The study of fluid queues has been the subject of much recent work. In these models the 
queue length is considered a deterministic (or "fluid" ) process, rather than a discrete random 
process that measures the number of customers. These models tend to be somewhat easier 
to analyze, as they allow for less randomness than more traditional queueing models. Also, 
the rougher description of a queue as a fluid is thought to be adequate for many important 
modern applications, such as ATM (asynchronous transfer mode) and other high speed 
integrated networks. 

Some of the earliest studies of fluid queues are due to Kosten [21] , (23 and Anick, Mitra 
and Sondhi |2j. We briefly describe the model in 0, since much of the latter work can be 

*ddominl@uic.edu 
tknesslQuicedu 



1 



viewed as extensions/generalizations of it. There are identical, independent sources that 
turn on and off at exponential waiting times. The rate from off to on is A, and time is 
scaled so that the rate of the reverse transition equals 1. When a particular source is on, it 
generates fluid at unit rate to a buffer. We denote the buffer content at time t by X{t). This 
buffer is depleted at a constant (deterministic) rate c, as long as it is non-empty. 

If c < A^, which we assume henceforth, the buffer content may be non-empty. Letting 
Z{t) be the number of on sources at time t, the buffer evolves according to the law 

• ^ / Z{t) - c, if X{t) > 

\ max [Z{t) - c, 0] , if X{t) = 

where [cj denotes the integer part of c. For simplicity we assume that c is not an integer, 
so that the fractional part of c, {c} = c — [cJ G (0, 1). 

Since there is randomness in the evolution of Z{t), the process {X{t),Z{t)) is stochas- 
tic and the exponential distributions of the on and off periods imply that the process is 
Markovian. These models are called Markov-modulated fluid flows. Of primary interest in 
applications is the distribution of the buffer content X{t). This can be obtained easily if we 
know the joint distribution Fk{x,t) 

Fk{x,t) = Pr [X{t) < X, Z{t) = k I X(0) = Xq, Z{0) = ko] 
0<x,xq; 0<k,ko<N 

but the latter is often difficult to compute. It seems that it is known explicitly only for the 
simplest models, such as those in |24j, [23] and P). Surveys of (single buffer) fluid models 
appear in [211, 03 and |29|. 

In j2] the steady state limit t ^ oo is considered for the model described above and the 
authors obtained an explicit expression for Fk{x) = Fk{x,oo), as a spectral representation. 

Then, the marginal buffer distribution was obtained from 2_^f^ Q^k{x). The complexity of 
this function has led to various asymptotic investigations. These assume that the number 
of sources N — > oo, with an appropriate scaling of the state variables {x, k). In [33] and jSH] 
the theory of large deviations is used to construct an approximation to the probability that 
the buffer exceeds the value x = Ny, of the form Pr [X(oo) > Ny] exp [— I(l/)] ? where 
l{y) satisfies a variational problem that can be solved fairly explicitly. 

A more complete asymptotic description of this buffer overflow probability is obtained in 
[HT] . where the author shows that 

Pr [X(oo) > Ny] ~ N^hiiy) exp [-N lo{y)] 

and it is also shown that lo{y) is equivalent to the solution of the variational problem in |35j. 
It is furthermore established that the asymptotic approximation is quite accurate numeri- 
cally. The technique used is to expand the spectral representation in |2I using asymptotic 
methods, such as the Euler-McLaurin formula and Laplace's method 

The existence of a steady state for this model requires that < 7 < 1, c = A^7, 
which says simply that the processing rate exceeds the average input rate. The case where 
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this stability condition is only weakly satisfied is called "heavy traffic". More precisely, this 
corresponds to the scaling 7 = + O (A^"^/^) . In j22I the authors obtain an approximate 
diffusion model in the heavy traffic limit and solve it explicitly as a spectral representation 
involving Hermite polynomials. We note that the discrete model has Fk{0) 7^ for < A; < 
\_c\ , which says that there is a non-zero probability that the buffer is empty for this range 
of on sources. Then Fk{0) = for k > c and this boundary condition was used to obtain 
explicitly the coefficients in the spectral representation in 

In the heavy traffic limit the joint distribution satisfies a parabolic PDE, that behaves as 
the heat equation in a part of the domain and the backward heat equation in the remainder 
pn] , [22] • Such problems arise in a variety of applications, such as counter-current separators 
[TB] , mean exit times ^21, the Milne problem of statistical physics jB], neutron transport 
theory J7] and diffusion in spatially varying convection fields. 

Their study goes back to Gevrey |T3I, [Hj and more recent analyses appear in [3], 
0; and |19j. The interesting mathematical feature of these problems is that the initial 
(or boundary) conditions can be imposed only where the PDE is forward parabolic. This 
"half-boundary condition" makes the problem difficult to analyze. The model in j2| may be 
viewed as a discrete analog of these PDEs and the "half BC" corresponds to the condition 
Fk{0) = 0, that can be applied only for [cj + 1 < A; < A^. 

In 120] we developed an asymptotic approach to analyze backward-forward parabolic 
PDEs, in a limit where the diffusion coefficient is small. It is based on the ray method of 
geometrical optics |TH] and matched asymptotic expansions. An important feature is the 
careful treatment of the point on the boundary where the PDE changes type from backward 
to forward parabolic. In this "corner region" the asymptotic solution may be represented as 
a contour integral involving Airy functions. 

The purpose of this paper is to extend the asymptotic approach in ^| to discrete models. 
We shall analyze the model in j2] directly by using the differential- difference equation satisfied 
by Fk{x). We make no recourse to the exact spectral representation of the solution given in 
[2]. After appropriate scalings of k and x, we shall analyze this equation asymptotically for 

EN 
Fk{x) in 
fc=0 

pi] can be easily recovered from our two-dimensional results. There are several important 
differences between the asymptotic structure of these discrete models and the backward- 
forward parabolic PDE studied in [20]. For example, the structure of the solution in the 
corner region, where a; ~ and k k, is different than the corresponding range in the 
diffusion model. Here the solution can be expressed in terms of Bessel functions. 

A variety of extensions/generalizations of the model in |2j have been recently analyzed. 
For example, the transient solution (more precisely, its Laplace transform over time) is 
obtained in [22 as a spectral expansion, using arguments similar to those in j2]. In [2Zj the 
author allows the fluid input rates from the on sources to depend upon the buffer size. This 
allows for an admission control policy for the fluid level (buffer size). An asymptotic analysis 
was done in which assumed that A^ is fixed, but that the input rates vary "weakly" with 
the buffer size x. Finite buffer size models are considered in j2H] and in [20] the model allows 
for certain sources to increase the buffer while others deplete it. The latter leads to a three 
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dimensional problem, where one must keep track of the buffer content and also the number 
of active sources of either type. Problems with two buffers and various priority mechanisms 
are studied in [Hj, [H] and |2l]. In jH] and ^H] fluid models with more general birth-death 
modulating processes were analyzed. 

It seems that for the more general models considered |H], [0], [IHI, [H], (HDj, the 
solutions are not particularly explicit. They involve either solving systems of equations or 
computing the eigenvalues of matrices numerically. The merit of our asymptotic approach is 
that it yields relatively simple formulas, which are both easy to numerically evaluate and also 
provide numerically accurate approximations to the performance measures, even for moder- 
ate values of the large parameter A^. We use similar scalings as the large deviations studies 
j2Hl, [Slj, jnS|, but in contrast to these studies we provide the full asymptotic approximation 
and not just the exponential growth/decay rate (in A^) of the performance measure. 

We also carefully treat various boundary and corner regions of the state space 



and indeed we show that their analysis is needed in order to obtain the asymptotic expansions 
away from the boundaries. We obtain detailed results for the model in [2] and develop the 
methodology to treat other models of this type. 

The paper is organized as follows. In Section 2 we state the basic equations. In Sections 
3-7 we analyze these in various ranges of the state space (0). In Section 8 we recover 
the one-dimensional results in [HT]. Finally, in Section 9 we summarize and interpret the 
asymptotics. 

2 Problem statement 

In the model proposed by Amick, Mitra and Sondhi j2], a data-handling switch receives 
messages from A^ mutually independent information sources, which independently and asyn- 
chronously alternate between the on and off state. The number of on sources forms a 
birth-death process Z{t) with birth rate = A(A^ — k) and death rate fik = k, where the 
rates are conditioned on Z(t) = k. Each source is on, on average, of the total time. An 
on source transmits at the uniform rate of 1 unit of information per unit of time. 

The switch has infinite capacity, and stores or buffers the incoming information that is in 
excess of the maximum transmission rate c of the output channel. The drift = k — c gives 
the rate of increase of X{t) (the buffer content at time t) when the birth-death process is in 
state k. That is, the rate of change of X{t) at time t is rz{t), provided rz(t) > or rz(t) < 
and X{t) > 0. If the buffer has emptied at time t, it stays empty as long as the drift remains 
negative. 

Following [SS] we define 



{{x,k):x>0, 0<k<N} 
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The stationary probabilities pk of the birth-death process can then be represented as 



Pk 



1 



> TT 
^.7=0 ■ 



In order that a stationary distribution for X{t) exists, the mean drift o^^^^ should be 
negative 

A 



TV 



(A+1) 

which gives the stability condition 



A + 1 



N -c 



> TT,- 
^.7=0 ^ 



r, < 



Setting 
and 



A ^ c 



Fk{x, t) = Pr [X{t) < X, Z{t) =k]] t, a; > 0, < A; < 



(2) 



= 0, [0, iV], 

the Kolmogorov forward equations for the Markov process {X{t), ^if)) are given by 

+ ''^k-Q^ = Afc-i-ffc-i + Affc+i-^fe+i — (Afc + /Xfc) Ffe. 

For the stationary distribution Fk{x) = Fk{x, oo) with the above rates and drift, we have 
BF 

(k - c)-^ ^X[N-(k- 1)] Fk-i + (k + l) Fk+i - [X(N -k) + k]Fk, < k < N. (3) 

Moreover, if the number of on sources k exceeds c, then the buffer content increases and the 
buffer can't be empty. Hence, 



Fk{0) =0, [cj + 1 < A; < A^. 



Also, 



Fk{oo) 



N 



A^ 0<k<N, 



{1 + X)^\k^ 

since Fk{oo) is the probability that k out of N sources are on simultaneously. 



(4) 

(5) 



5 



3 The ray expansion 

To analyze the problem dSj)-© for large we introduce the scaled variables y,z,'y, with 

k = zN, c = 'yN, X = yN, 2, 7, ?/ = 0(1). 
We define the function G{y, z) and the small parameter e by 

e = ^, Fk{x) = G{ex,ek) =G{y,z) 

and note that Fk±i{x) = G{y, z ± e). 

Then Q becomes the following equation for G{y, z) 

e{z - i)Gyiy, z) = X{l-z- e)G{y, z - e) + {z + e)G{y, z + e) - [A(l - z) + z] G{y, z) (6) 

and (jlj) implies that 

G'(0,z) = 0, 7<z<l. (7) 

To find G{y, z) for e small, we shall use the ray method. Thus, we consider solutions 
which have the asymptotic form 



G{y, z) ~ £^ exp 



-^{y,z) 



^{y,z). 



Using © in ©, with 



U!{y, z±e) = U!±^, + Uf^^s + O {e^) 



dividing by exp [^^^(y, z)\ , and expanding in powers of e we get 



[Z - 7) {^<yK + 



+ 

where 



zU + \{z~l)- 



zU + {X-l)z ~ \ + \(l - z)- 



1 



1 + ) f/ + A ( 1 + - 1 ^ 



Uiy,z)=exp [^,iy,z)]. 
Equating the coefficients of e we obtain the eikonal equation for ^(y, z) 

{z - 7)^y + {1-X)z + X + X{z -l)^-zU = 



and the transport equation for K(?/, z) 



l + ;7^.Jf/ + A 1 + -$ 



1 

f7 



u 



K + (7 - ;z)Kj, 



z?7 + A(z- 1) 



(9) 



(10) 



0. 
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To solve (jlOj) and (fTT|) we use the method of characteristics, which we briefly review below. 
Given the first order partial differential equation 

d {y,z,^,p,q) = 0, 

where p = "^y, q = "^z, we search for a solution \E'(?/, 2;). The technique is to solve the 
system of "characteristic equations" given by 

y — -g^ — Vp, Z — 

p = -dy - pd^, q = -dz - qd^ 
^ = pdp + qdq 

where we now consider {y, z,ipjP,q} to all be functions of the variable t, with ip{s,t) = 
^{y,z). 

For the eikonal equation (jiup . the characteristic equations are 



y = z-j (12a) 

i = A(l - 2)6-" - ze'^ (12b) 

p = (12c) 

q = e'i - Xe-1 -p+ \-l (12d) 

^ = p(z--f) + q[\{l-z)e-''-ze'^]. (12e) 



The particular solution is determined by the initial conditions at t = 0. We shall show that 
for this problem two different types of solutions are needed; these correspond to two distinct 
families of characteristic curves, or rays. 

3.1 The partial derivatives "^y, 
Setting "^ylf-^Q = s, and solving ()12c|) yields 

p = s, 

so that is constant along a ray. Introducing the function u{s,t) as in Q 

u{s,t) = exp[q{s,t)] 

we have from ()12d|) 

f ^ . \ 

ii = uq = { u s + A — 1 ] u 

\ u J 
= u^ + s)u-\ (13) 

= [u- ri){u - 
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where 

1 



riM = -{s + l-X±A), A{s) = ^{X-s-iy + 4X 



and ri corresponds to the (+) sign. 
Solving (fT^ gives 



{u - ri) (mq - r2) 
{u - r2) {uq - ri) 



t, u{s,0) = Uo{s) (14) 



or ^ 

u{s, t) = r2 H 7 r . (15) 

gAt ( ^ - 1 ) + 1 

Evaluating (fTUI) at t = we get 

(1 - A)7 + A + A(7 - 1)— -juo = 

Uq 

so that ^ 

Mo = 1 or Mo = - (1 - 7) . (16) 

7 

3.2 The rays from (0, 7) 

We now consider the family of rays emanating from the point y = 0, z = 7. From (jl2bj) 
the equation for z is 

z = X(l~z)--zu (17) 
u 



which we can rewrite as 



whose solution is 



dz . 1 

—u = Xil — z) zu 

du u 

C{s)u-X 



u^ + {X-l-s)u-X' 

Imposing the initial condition z(s, 0) = 7 and using ()16|) . we obtain C{s) = A — 7s for both 
possible values of uq. Therefore, 

(A — 7s)?i — A 



u^ + {X-l-s)u-X' ^ ' 

From p2|) the equation for y is 

y = z-'y (19) 

which implies that 

yis,0)=zis, 0)-7 = 



8 



and 



y{s, 0) = z{s, 0) = A(l - 7)— - 7Mo- 

Uo 



We define 

p = 7 - A + A7 
with < p < 1 from (j2I). From (fT^ we have 



(20) 



y{s,o) 



-p, Uo = l 

p, uo = H^-l) ' 



Using the initial condition y{s, 0) = and expanding in powers of t, we get 

y{s,t)r~.y{s,0)-, t^O 

and in order to have y > for t > (i.e., for the rays to enter the domain [0, 00) x [0, 1]) 
we need to choose 

«o = -(l-7) (21) 

7 

with Mo < 1 from 

Integrating (fT^ and using (fT^ . (fTij) and (fTH)) . we have 



y ~ J -^(■^5 '^^)c^^ ~ 7^ 


u 

(A — 7s)w — A 



uo 



[w2 + (A - 1 - s) w - A] 



— 7— In 



{u - ri) {uq - r2) 
{u - ra) (mo - ri)_ 



-i^[(A + l)p + <^.]ln 



(M-ri)(Mo-r2) " 

(M-r2)(Mo-n). 

[r2(A - 57) - \]{u- Uo) 



where 



^ [ri(A - 57) - A] (u-uq) ^ 

A2 (m - n) (mo - ri) A2 (u - r2) (mq - ^2) 

= 7 + A — 7A 



with > from 

Finally, combining (US)), (UHl), (ES) and (ED) we conclude that 



yis,t) 



1 

A2 



bs + {X + l)p . 



(22) 



(23) 



A 



sinh(At) + pcosh(At) - p + [75^ + (7 - A - A7)s + A(A + 1)] t| 

(24) 
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Figure 1: A sketch of the rays from (0,7). 



z{s,t) = — [0s + (A + l)pcosh(At)+pAsmh(At)+7s^ + (7-A- A7)s + A(A + 1)] . 

(25) 

This yields the rays that emanate from (0, 7) in parametric form. Several rays are sketched 
in Figure 1. 

Combining the terms depending on t we have another useful expression for y{s, t) and 
z{s,t) 

yi^^ t)-^ (^^* - 1) + ^ (1 - --^') + [^3(.) - 7] t (26) 

z{s, t) = A,{s)e^' + A2(s)e-^* + A^is) (27) 

with 

A,{s)^^{X + l-s + A) (S7 + P-A + 7A) 

Ms) = ^ [A(l + A) + (7 - 7A - A) s + 7s'] . 
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For t > and each value of s, (j^^ and ()25p determine a ray in the {y, z) plane, which 
starts from (0,7) at t = 0. We discuss two particular rays which can be obtained in an 
explicit form. For s = we can eliminate t from p5|) and obtain 



1 



A + 1 



In 



z\ + z — \ 
P 



7 < z < 1 



(28) 



and using we get 



y = Yo{z) 



z — ^ 



A + 1 (A + 1)^ 



In 



z\ + z — \ 
P 



s = 0, 7 < 2 < 1. 



(29) 



Also, from ()18|) we have 



\ — z 

u = A , for s = 0, J < z < 1. 



(30) 



For s > Sq = ) we have both y(t) and z(t) increasing for t > 0. When 

s = Sq the function u is constant along the ray 

u{So,t) = Mo = - (1 - 7) , 

7 

2 

y(t) increases and z{t) asymptotes ^ with 

5= (7-1)2a + 72, S>0. 

Eliminating t we obtain 

"7(1 -7)p 



t = ^^i^ln 



7^ — 



7 



7<.<y 



y = ri(^) 



[7(1 -7)] V 
52 



7^ — 5^; 
.7(1 -7)P 



+ In 



" 7(1 -7)p " 
7^ — ^2; 



- 1 



s = So, 7 < z < 



7 



(31) 



For s < Sq the rays reach a maximum value in z at t = Tmax, where 



J- rr 



and we have 



|/(s,Tmax) = --r^ + 



2 A 1^1 



P , ^3 - 7 



A2 ■ 2A 



(32) 
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Figure 2: A sketch of the rays Yq{z),Yi{z) and the curve Fmax(-2^)- 



Inverting ()33|) we find that at t = Tmax 

_z + X{l-z)± 2^zX{l - z) 

7 — 2; 



(33) 



(34) 



From (jH^ we see that ^ — as s t Sq. Therefore we must choose the (— ) sign in (jH^ 
and conclude that 



7 — 2; b 



which when used in yields 



y ^^ax(2^) 



p 



+ 



^3(^n 



7 



A2 (^max) 2A (5, 



In 



'^2 (5n 



(5n 



7 < 2; < 



r 
5 



(35) 



(36) 



This is the value of y when z reaches its maximum value. In Figure 2 we sketch the rays 
Yo(-2), Y\{^z\ the critical curve l^max(^) and a pair of rays showing their intersection with 
^max(-2) at the maximum. 
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Solving for t in ()27|) we obtain 

T+{s,z), <y < Y^s.x{z), 7 < 2: < 1 



T_{s,z), \y> Y^^^{z), ^ < z <^-^\u{y > 0, < z < 



(37) 



with 



T±(s, z) = —\n 



1 

A 



z-A3±J{z-A3Y-4A^A2 



2Ai 



Solving for u in we find that 



u 



f/_(s, z), 0<y < Fmax(^), 7 < ^ < 1 

f/+(s,2;), I?/ > Fmaxl^;), 7 < 2; < ^1 U {y > 0, < z < 7} 



(38) 



where 



U±is,z) = - s + l-A + -(A-7s) 
2 L z 

± ^ Vp' + [2 (A + 1) p + 20s] (2 - 7) + [(A + 1)2 + 2(1 - A)s + s^] {z - 7)2. 



From (|Tn|l we see that the maximum value in y is achieved at the same time that z = 'j, 



and that occurs at t = with 



1 



y ^ ^ >g-pA + (A+l)p ' 
^ A ^ [0s + pA + (A + l)p 

When -2 = 7, y > we also have from (|18|) 

m(s,T^) = 1, s < S'o, 

so that = when ^ = 7. 

Inverting the equations we can write 



s < So. 



(39) 



(40) 



s = S{y,z), t = T{y,z) 



and 



^{y, z) = tlj[S {y, z) , T(y, z)] , K{y, z) = K[S (y, ^) , T(y, ;2)] . 
We will use this notation in the rest of the article. 



13 



3.3 The function ^ 

From (I12ell we have 



ip = sy + \n{u)z 



which we can integrate to get 



dz 

iIj^s, t) = ip{s, 0) + sy{s, t) + / in [u{s, r)] — (s, r)dr 



iIjo^s) + sy + \n(u)z — In (mq) 7 — / —du. 



u 



Using (fTHj) and (fT^ we can write 



t 



—du 

u 



(A — 7s)u — A 
u(m — ri){u — 



du 



1 ^ (A — ■ys)ri — A (A — ■ys)r2 — A 



u ri(u — ri)A r2(u — r2)A 

" 

u\ (A — 7s)ri — A, f u — ri\ (A — 7s)r2 — A , ( u — 



du 



In 



Mo 



riA 



Mo - ri 



In 



Mo - ^2 



Hence, 



t) = V'o(s) + + (z - 1) ln(n) + (1 - 7) In (%) + ^^^."'^ ^ In ^ ""^ 



raA 



% - r2 



(A — 7s)ri ~ J ( u — r\ 



riA 



Mo - ri 



(41) 



Obviously, ■{/'(s, 0) = ■i/'o is a constant since all rays start at the same point. We will determine 
'?/'o in section 3. 



3.4 The transport equation 

Re-writing the transport equation ()11|) as 



l + |v,,.)f/ + A(l + lv,,.-|v,Jl 



lL={z- 7)K^ + 



zV ^\{z-\) 



U 



and taking into consideration (fT^ and (fTTj). we see that 



1 + 2^^^ j M + A I 1 + -i^zz - -Vzz 1 - 



K = yKy + zK, = K 
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or 



K 
K 



1 + 7;^V>zz ]u+- 
2 u 



1 + 



A 1 
u 2 



A 



zu-\ — (1 — 2;) 
u 



4>z 



\ ( AA 1 a 



zu (1 — z\ 

u 



since u = exp [^] . But, zu — ^(1 — z) = —z, so that 



K I A 
— ^ - [u + - 
K 2\ u 



Idz 

2dz 



l(u+-^ -- (z— + ——\ 



2 \ ' u J 2 \ dz dsdzj 
Introducing the Jacobian of the transformation from Cartesian (y, z) to ray (s, t) coordinates 

dz dy dz dy 



we have 



and 



Hence, 



ds 
dz 



dt ds ds dt 

y 1 - z 



(42) 



dt 1 dy 
dz J ds 



..dy ^ . dy dz . dz .. .. dy dz , . 
^ ds ^ ds ds^ ds^ ^ ds ds ^ ^ 



K I f A 

K = -2V'^u 

1 / A 
= - h^, + - 

2\ u 



1 


' ..dy 


dz , 


-7) 


2J 


'dS 


j 


1 


' A z 


j 


2J 


~ 2 


uz z 


~ J 



UJ z 
ijj z 



where a;(s,t) satisfies the ODE 
which can be solved to give 



d_ 
di 



ln(a;) 



uz 



U!{s, t) 



(A — 57) u — \ 
m[p + 7(1 -7)5]' 



(43) 
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We conclude that 



A-(M)^A.M.-y^, (44) 

with Kq{s) to be determined. As {y, z) (0,7) (i.e., as t 0) the Jacobian J(s,t) 0. 
Therefore, the asymptotic expansion in (jH)) ceases to be vahd. 

So far we have determined the exponent ip{s,t) and the leading amplitude K{s,t) 
except for the constant ipo in PT|l . the function Ko{s) in and the power z/ in (jHl). 
In section 4 we will determine them by matching (jHJ to a corner layer solution valid in a 
neighborhood of the point (0,7). 

3.5 The rays from infinity 

From (0), we have 

Ffc(oo) = ^(oo, z) = exp 1^ [z\n{X) - ln(A + 



and by Stirling's formula 

"1 



G(oo, z) ~ y/eK{z) exp 



^<l>{z) 



, < z < 1 (45) 



where 



^{z) = -z \n{z) - (1 - ^) ln(l - z) + z ln(A) - ln(A + 1) (46) 



Note that $(2) and k(z) satisfy (fTUI) and (fTT|) . respectively. 
Denoting the scaled domain by 

2) = [0,00) X [0,1], (48) 

we must determine what part of 2D the rays from infinity fill. Defining p^o = ^ and 

Moo = exp (H ) , we have 

= 0, Moo = ^ -■ (49) 

z 

(fT^ and ()12b|) yield the equations for the rays yoo(t), Zoo(t) 

yoc = Zoo- 7, ioo = (1 + A) 2:00 - A (50) 
or, eliminating t from the system and writing yoo{t) = Yoo{z) we get 

dYoo z-j 



dz {1 + X)z-X' 
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(51) 




Figure 3: A sketch of the rays from infinity. 
Solving subject to the initial condition Yoo{zq) = i/q, where (yo, Zq) G dD we get 



Y^{z) = yo + 

A 



z- Zo 



In 



(A + l)z- A 
(A + l)zo - A 



A + 1 (A + r' 
and when zq = the ray is a line parallel to the y-axis given by 



Zo ^ 



A 



A + 1 



(2/00, Zoo) = [y 



A 



A + 1 



0<y. 



(52) 



(53) 



We can interpret ()52|) as a family of rays emanating from the "point" (00, ^^) and hitting 
the domain's boundary dD at the point {yQ,zo). We divide the rays into five groups (see 
Figure 3) depending on the location of (?/o,-Zo)- The quantities inside the logarithms are all 
greater than one: 

1) yo > 0, zo = 0. 

A 



Y^.(z) 



yo + 



+ 



p 



A + 1 ' (A + 1)' 
2) yo = 0, < Zo < XTT- 



In 



A - (A + 1)0 



Yoo(z) 



z - Zo 



+ 



A + 1 (A + 1)^ 



In 



(A + 1)^0 - A 
_{X + l)z-X_ 
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0< z < 



Zq < Z < 



A 

aTT' 



A + 1 



For zq = the ray becomes parallel to the y-axis 
3) 2/0 = 0, ^ < 2;o < 7. 



Z- Zq 



P 



A + 1 (A+1)^ 



In 



(A + l)^o- A l A 
_(A + 1)z-aJ' A + 1 



< Z < Zq. 



4) ?/o = 0, ^0 = 7- 



z — 7 



A + 1 (A + 1)^ 



In 



P 



{X + l)z-X 



A + 1 



< 2 < 1. 



This critical ray is tangent to the z-axis at z = 7. The upper branch of the ray for 
7 < 2; < 1 is the same as the ray Yq{z) which emanated from (0, 7) and was defined in ()29|) . 

The part of the boundary corresponding to {0} x (7, 1] U [0, Yo{l)) x {1} is not reached 
by the rays. 

5) vo > Yo{l), zo = 1. 



Yooiz) 



yo 



p 



In 



(A + l);z-A 



A 



A + 1 



A + 1 (A + 1)' 
The rays from infinity fill the region given by 

R = {0<y, < 2 < 7} U {Yo{z) <y, 7 < ^ < 1} 
The complementary region R'-^ 

R^ = {0<y<Yo{z), 7<^<1}, 



< z<l. 



(54) 



(55) 



is a shadow of the rays from infinity. In i?*", G is given by (jH)) as only the rays from 
(0,7) are present. In the region R below, both the rays coming from (0,7) and the rays 
coming from infinity must be taken into account. We add (jHl) and to represent G in 
the asymptotic form 



G{y, z) ~ v^k(z) exp 









U{z) 


+ exp 











^{y,z), {y,z)eR. 



(56) 



We can show that $(-2) > "^{y, z) in the interior of R, so that G{y, z) ~ G{oo, z). However, 
in R we can refine to G{y, z) — G{oo, z) ~ e'^ exp [^^(|/, z)\ K(?/, z). 



4 The corner layer at (0, 7) 

We determine the constant z/ in (jH)) and the function Kq{s) in (j44j) by considering carefully 
the region where the rays from (0, 7) enter the domain 25 and using asymptotic matching. 
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We introduce the stretched variables x, 1, the function Gi{l, x) and the parameter a defined 
by 



X = p X>0, Fk{x) = Gi{k-c + a,^) = Gi{l,x) 

I = k — c + a, — oo < / < oo 
a = c — [c\ , < « < 1. 

Note that a is the fractional part of c. Use of (jHTjl in (jH)) yields the equation 

(/ - «)^^"' = bGi + 1, X) + A(l - 7)^1 (/ - 1, x) - 0^1 (/, x)] e-' 
+\{a + 1 - l)Gi (/ - 1, x) + (/ - a + (/ + 1, x) + (A - 1) (/ - a)Gi (/, x) 
or, to leading order in e, 

^(1) 



(57) 



(/-a)- 



7F£1 + A(l - i)Fl^[ - ct^F, 



(1) 



dx 



(5^ 



with Gi (/, x) ~ i^/^^ (x) as £ ^ and = 7 + A - 7A. 
Also (jH) gives the boundary condition 

f/^)(o) = o, />1 

and (gil) implies that Ffc(oo) ~ fI^\oo) 

fI^\oo) = {uof^" ^(7) exp 



i$(7) 



(59) 



(60) 



with Mo = ^ (1 — 7) . For a fixed / and x ^ 00, we approach the interior of R, where (j^Bjl 

applies. Note that f/"'^'*(oo) is the same as ^/eK{z) exp [^<l>(z)] expanded for z 7. Thus, 
(jUn|) is the asymptotic matching condition between the corner layer and the solution in R. 
We shall examine the matching to i?*" later. 
Equation (j58|) admits the separable solutions 



if hi{6) satisfies the difference equation 

7/ii+i + A(l - 7)/ii_i = [(/ - a)e + 0] hi. 
Setting hi{9) = (uq)'^ Hi{6) we see that 



Hi+^ + Hi_^ = ^ [(^ - + ^] Hi 



(61) 



(62) 
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with 



/3 = 2v/A7(l-7)- (63) 
The only solutions to ()fi2|) which have acceptable behavior as / ^ oo are of the form 



where J is the Bessel function. If is not to grow as x — > oo, we need < 0. But except 
when V is an integer, the Bessel function Jt,(x) is complex for negative argument. Therefore, 
we need 

^ - a = -1,-2,... 



9 



or 



<0, j>0. 



' J + 1 - a 
It follows that the general solution to (j58p takes the form 



(64) 



j>0 



or 



f[^\x) = v^K)' "fi:(7)exp 



l-a 



+ iV^f exp 
i>o 



-Hi) 

e 



(65) 



j + l-a 



X Ji-i-j 



-?(j + l-«) 



where the coefficients Qj in the above (spectral) representation remain to be determined. 
Taking the Laplace transform 



of we obtain 



F«(^) = v^(«o)'"°«:(7)exp 

1 







i$(7) 



1 



j>0 *^ + i+l-a 



t(j + i-«) 



(66) 



Thus, the only singularities of f/^''(^9) are simple poles at -i? = and = 6j, j > 0. It 
is well known that the gamma function r(2;) has simple poles at z = 0, —1, —2, . . .. Hence, 
we shall represent Fi{'&) as 



/3 



(67) 
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where /("i?) is chosen such that F (| + 1 — a) Ji_^_^_i (^) /("i?) is analytic for Re('t9) > - 
Taking the Laplace transform in we get the equation 

{I - a)^F^^^ = ^F^l\ + A(l - 7)f/_'} - / > 1 

which is satisfied by fl67|) . By the inversion formula of the Laplace transform we have 



27ii 



Br 



where Br is a vertical contour in the complex "i?— plane on which Re('(9) > 0. 
Since the residue of F^^\i}) at i) = corresponds to F^^\oo), we must have 



as -i? — > 0. In that limit we see that 



r(5 + i-«)>/^ 



Therefore, we write 



f{'d) = Ve{^) "./■^K(7)exp 



-$(7) + T(^9) 



where 



T{d) = ( ^ - a 1 In 



^\ , 2A(l-7) , . ^ 



and f{'d) is entire, with /(O) = 1. 

By combining the preceding results we have 



Fi'\x) = V~eJ^n{l){V^ot''e^^ 



IP 



-$(7) 

£ 



2m 



^"41 ? + J, 



Br 



The boundary condition (|^ implies that 



lim 

1?— >oo 



^fI^\^) 



/3 



= 0, l>l. 



exp [T(?9)] 



From (|7U|) we have 
1, 



+ 1 - a J; 



exp [T(?9)] 



F(l-a) ((3 



d) F(/-« + l)V2 



l-a 



0", ^ 
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Fixing / = 1, we get f{'d) = o {"d) , ^9 — » oo and by Liouville's theorem this forces f{^) to 
be a constant. Since /(O) = 1, we have f{{}) = 1. Thus, flTUj) becomes 



l-a 



-Ml) 



(71) 



2m J ^ \^ 

Br 



exp [T{^)] (I'd. 



The coefficients aj in the spectral expansion ()65|) are determined from ()71|) by applying the 
residue theorem. Noting that 



Res 



r(^ + i-a), 



we obtain 



aj = y/e. —K{'y) (a/uq) "exp 



-$(7) + Tie, 



U + l-aY J\ 



j\ (j + 1 - a) 



J > 0. (72) 



This completes the determination of the spectral and integral representations of f/^"* (x) and 
hence the leading term for Fk{x) in the corner region. 



4.1 Matching the corner and regions 

In this section we shall determine the constant ipo in ()41|). the function Kq[s) in and 
the power u in (jSJ. From we have 



If we introduce the new variable Q defined by 



n 



2py 2px 



(2: -7)2 {I -ay 



Q = 0(1) 



(73) 



we see that the limit Xy I ^ 00, Q > 1, corresponds to the matching between the corner 
and R'" regions. 
We set 

^ = ee, r] = {I - a)^ + (f) = {z ~ + (f), 77,6 = 0(1), r/, 6 > 
use this in Ji_^^± (^) and let e — > 0, thus obtaining 



exp < V 



J I - {(3/v f + In ( 



2Tiv^i-{PMy 
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with 



t; = p{r]) = ^7]^- (3^, p{(j)) = p. 



Use of Stirling's formula gives 



and from ()69|) we have 

exp [T(e9)] = exp 
Therefore, 



: exp 



1 

7e 



V 

1 

7q 



>\n (^^j +2A(l-7)-^lnM 



£0 



■^M ) r ( ^ + 1 - « ) exp [T{eQ)] 



P{v) — In 



Using (f73j) in (f7T|) yields, in terms of z and f2 



f[^\x) ~ v^yp«;(7) exp <j ^ 



/5 ; 



p - - In (uo) 



X 



27ri 



1 



Br' 



exp 



dr] 



where 



9{r]) = 7' + 



P J 



p - - In (uo) 



(74) 



(75) 



(76) 



(77) 



2p ^ ~ 

and Br' is a vertical contour in the complex plane with Ke{ri) > (p. For e —>■ with Q fixed 
we can evaluate (fTUj) by the saddle point method to get 



Fi^\x) ~ v^VPfi:(7)exp 
e 1 1 1 



\ $(7) + ^Y^ln(no) 



exp -(2; -7)^(77*) 



(7J 



'g"{ri* 



where the saddle point t]* (Q) is defined by g'ijf) = 0. Note that rj* (fl) > for f2 > 1, i.e., 
the saddle point 77* (f2) lies to the right of the pole at 77 = and the integrand is analytic 
for Re (77) > 0. 
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Taking derivatives in (f77j) we find that 



9{v) = 7r + 



p + 



9"{V) 



2p {ri 



r/2 + 2r](j) - 02 _ 2/52 
p{r]) 



|lnM + 01nr^^Vp(^) 



\ P J 



- 20 In 



V /5 ; 



2p - 01n (uq) 



(79) 



(80) 



From ()79p we observe that g'{i]*) =0 if rj* = cj) and = 1, which implies that 
Tj* (1) = 0. To determine 77* for ~ 1, we use fl79p and an expansion of the form 



T]* (fi) ~ + ai (fi - 1) + 02 (^^ -if + a^{VL-lf + ■■■ . 
Using (jHT|) in (f?^ and expanding the latter in powers of — 1, we find that 



^1) 



and 



^^(r;*) ~ 1 In (no) - g (r^ - 1)\ g"{v* 



(fi-l)-\ 



1 



r]* - 3p2 ' ^p(ry*) ^ 

from which we conclude that as (y, z) —>■ (0,7) 

F/^)(v) ^ _ Ai L 



L= (1] - 1)-^ exp I i [$(7) + (^ - 7) In («o)] } . (82) 



We next evaluate K and -0 in (jH)) near the corner (0,7). Using T+{s,z), (cf. dSZI)) in (EH) 
and expanding for small s, we obtain 



with 



y^Yo{z)-Y2{z)s 
2C 



z — 7 



(A+1) (A;2 + z-A)^ 



(A + 1)^ 
2Cp 



In 



z + zX — \ 



P 



(83) 



M) 



+ 



3C 



(A + 1)^ (A + 1) 



(;.-7) + (A-l)(^-7)' 



and 



C = 2A-7+(7-l)A2. 



^5) 



When z is close to 7, (jH!^ yields 

y 



1 / x2 (A + l)p + 0s 3 
^(^-7) p (^-7) , ^-7 
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or, using 



2 2; — 7 



(A + 1)^, 



Also, from ()37|) we have 



2; — 7 
t ~ , 2 — > 7. 



We expand PT|) for small t 

^p{s,t) r-^ ^p(j + \■a{uo) pt, t^O 
which taking (|H7|) into account gives 

'^{y,z) ipo + (z - 'j) In (uo) , 

in agreement with ()82p if 

^0 = $(7)- 



From we obtain 



K{s,t) K,{s) — 



t 



or, using (j86|l and (jSTjl in the above. 



K{y, z) ~ J^o 



2 z — 7 



P 



7 1/(1-7)0(^-7) 



3 ' 



2; ^ 7. 



Matching the algebraic factors in (jH2I) and (jEHI) yields 



which implies that 



3p2fi- 1 
2 (f) z — 'J 



Ms) 



1 (f) Z — 'J 

3n pi Q — 1 



27rs' 



(86) 
(87) 



(89) 



(90) 



Finally, the exponent of £ in (jH)) is determined from (j82|) and turns out to be u = 1. 
This completes the determination of the asymptotic solution corresponding to rays from the 
point (0,7). To summarize, we have established the following. 



Result 1 The solution of (0j is asymptotically given by 

'1 



G{y,z) ~ eexp 



:'^iy,z) 



]K(?/, z) in R 



c 



(91) 
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with 



G{oo, z) — G{y, z) ~ — e exp 



K(?/, z) in R 



G(cx), z) = exp <j i [2 ln(A) - ln(A + 1)] | f 



= {0 < < 2; < 7} U {ro(^) < 7 < ^ < 1} 



Y.iz) 



z — 7 



In 



A + 1 (A + 1)' L(A + 1)2;-A_ 



, 1 <z<l 



A 



p = 7A + 7-A, Mo = -(1-7) 

7 



(92) 



^'(y, z) =^(s,t) = ln 



A 



7 (A + 1), 



+ s?/ + (2;- l)ln(M) + 



(A - 7s)r2 - A ( u-r2 



In 



^^0 - r2 



(93) 



(A — 7s)ri — A f u — Ti 



riA 



In 



uo - ri 



(94) 



ri,2(s) = ^ (s + 1 - A ± A) , A(s) = v/(A-s-l)V4A, 



(?/, z) is related to {s,t) by |^j), ^EW and u{s,t), u!{s,t), 3{s,t) are defined in i fT^)) . (|73D 
and respectively. We note that s < in R so that the right side of KU^) is positive. This 
gives the leading term for the probability 



Pr 



y z 
X(oo) > X = -, Z(oo) = k = - 



that the buffer exceeds x = Ny. 

In the corner range where \57\ ) applies, the leading term is given by ^71]) , or ( t6'5)) and 

5 Transition layer 

We shall find a transition layer solution near the curve y = Yo{z) defined by ()29|1 which 
separates R and i?*" . On this curve s = 0, hence is not valid because K{y, z) is infinite 
there. 

We introduce the new function L^^x) defined by 



FJx) = Fk((X))Lk(x). 
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Then Q yields for Lk{x) the equation 

{k - c)L'„ = [X(k-N)-k]Lk + X{N- k) Lu+i + kL^-i 
and matching to region R forces 

Lk{^) = 1. (95) 

In terms of the variables y = ex, z = ek, the function L^^\y, z) = Lk{x) and the parame- 
ters 7 = ec, e = j^, we get 

.dim aid) £ a^-Lm ^fcfiim , 

Introducing the stretched variable A defined by 

y = Yo{z) + y^A (96) 

and the function L^'^\A, z) = L^^\y, z), we obtain for L^^)(A, z) to leading order the diffusion 
equation 

(,_,_,,)______,,_,)^_.0. (97) 

To solve (jnZj) we assume that L'^^)(A,z) is a function of the similarity variable V = -j^, 
and let ii(V) = L^'^\A,z), where /i(z) is not yet determined. From we get 

-i-' (A2; — 2; — A) (2; — 7) 

and gives 

£(00) = 1. (99) 
We can eliminate A by choosing fi to satisfy the equation 

o / (A^ + ^-A)^ 

2/i/i' ^ 2 = -1 

(Xz — z — A) (2; — 7) 

or equivalently 

d{fx^) {Xz-z-X){z--fy 



— - q . (100) 

dz ^Xz + z-Xf 

We choose /i(7) = 0, which is necessary for matching the transition layer with the corner 
layer solution (fTTj) . and solve pOO|) to obtain 



fj.{z) = (101) 

where Y2(z) was defined in 
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Now 



and are just 



£" = -VS.\ £(oo) = 1 



and the solution is 



V 



Z{V) 



Thus, the transition layer solution for y — Yq{z) = O (v^) and 7 < z < 1 is 

V 

Fk{x) ~ F^'^\V, z) = ^/eK{z) exp 



)= I exp ] dr. (102) 



5.1 Matching the transition layer and solutions 

We show that as ()102|) is expanded as V 



y-Yo{z) 



-00, the transition layer matches to 



the ray expansion, corresponding to rays emanating from (0,7). As ?/ Yo{z)-, from 
and (1^ we have 



y - Yoiz) 



ny,z] 



1 , fzX + z-\ 
In 



Y2{z) ' A + 1 V P 

Using (jl03j) in (jU^ and and expanding for y Yq{z), we obtain 

i[y-Yo{z)f 



^iy,z) ~ <l>(z) 



2 Yo(z) 



K{y,z) 

From ()102|) we have 

F(2)(F,z) ~ -v/^exp 
in agreement with (jl04p and (jlOSp . 



1 WY,(z] 



27ry-Yoiz) ^ z{\ - z) 



z\ 



exp — 



V 



-00 



(103) 

(104) 
(105) 



6 The boundary layers at z = and z = 1 
6.1 The boundary layer at z = 

From (j94|) we see that K{s,t) is singular as z 0. Therefore, we find a boundary layer 
correction near z = 0. We consider solutions of which have the asymptotic form 



Fk{x) - Ffc(oo) = Flf\y) - Ffc(oo) ~ ^^^-'^exp 



-^(y,0) 



K!^\y)- 



(106) 
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Using p06p in Q and expanding in powers of e gives to leading order 



Solving ()l()7j) we obtain 



Q = {k + l)Kfl + [^^y{y,Q)-X]Kf\ 



and hence 



VQ-k 



-^(2/,0) 



(107) 



(108) 



(109) 



Setting k = z/e, F'f'\y) — Fk{oo) = G^^\y, z) and letting k ^ oo, we get 



G^^\y,z) ~£'^Oexp|J*(?/,0) + ^zln 
From ()92j) we have 



A- 7^,(1/, 0) 



£ 1 



G{y, z) - G(oo, z) ~ ^ ^^^^^ exp |ivl/(^, o) + -z In 

27r V-2o(2/, 0) ye e 



27r 



A -7^(^,0) 



k^^\y). (110) 



A- 7^(1/, 0) 



[7(1- 7)^(1/, 0) + p] Jo(z/)' 



z ^ 



(111) 



with 

5(y, 0) = ^y{y, 0), Jo(2/) = bS{y, 0) - A] ^ [S(y, 0), T(y, 0)] + 7^^ [5(y, 0), T(y, 0)] . 
Matching (jllOj) and (|11H) we conclude that 



.. = i /,(3).,_ VP I A-7g(y,0) 
'2' ^s{y,0)y [l{l-l)S{y,0)+p]3o{yy 



Therefore, 



F^^\y)~Fkioo)^e^^-'exp 



-m{y, 0) 



[A - lS{y, Q)f ^ 



A- 7^(1/, 0) 



V2^S{y, 0) V [7(1 - l)S{y, 0) + p] Jo(z/) ' 
We note that the the right side of ()112|) is negative since Siy, 0) < < 0. 



(112) 
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6.2 The boundary layer at z = 1, < y < Yq{1) 
From (j9H) we see that 



VP 



-[l + (l-7)g(l/,l)] 
27r5(y,l) V [p + i {I - ^) S{y,l)]J,{y) 



{l-zY 



z 1 



(113) 



where 



UV) = [(1 - 7) Siy, 1) + 1] ^ [S{y, 1), T(y, 1)] - (1 - 7)|^ [S{y, 1), T(2/, 1)] , 

so that ]K(?/, 2;) is singular when z = 1. Therefore, we introduce the new variable j = N — k 
and consider solutions to (jHl) of the form 



Using ()114p in Q gives, as e ^ 0, 



(114) 



[1 + (1 - i)%iy, 1)] - A (j + 1) <\ = 



which we can solve to obtain 

'(4) 



l + (l-7)vl/^(y, 1) 



A 



Hence, 



From we have, as 2; — 1, 



A 



(115) 



^{y,z) ~ ^(i/,l) + ln 



1 + (1-7)5(|/,1) 



(l-z), Siy,l) = ^yiy,l)- (116) 



A(l-z) 

Expanding (|115|) as j — >^ 00 using Stirling's formula and j = yields 

■1 + (1- 7)^,(1/, 1) 



Fj'\y) ~ e--^ ^ J^ k^'\y) exp (-v&(l/, 1) + - In 



/27r VI - ^ 

Matching (ITTT^ and (fTTT)|l with (fTTTfl we conclude that 



A(l-2) 



(l-z) 



(117) 



1 ^^.(4)., VP I -[l + (l-7)g(l/,l)] 

'2' x/2^5(?/,l)V[P + 7(l-7)^(z/,l)]Ji(y)' 
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Therefore, 



1 + (1-7)^(1/, 1) 1^' 1 
A J J\ 



:ii8) 



X 



[l + (l-7)^(l/,l)] 



^S{y, 1) V [p + 7 (1 - 7) S{y, 1)] Ji(y) ' 
In the range < y < lo(l)) we have S'(?/, 1) > 0. 

6.3 The corner layer near {Yq{1)^ 1) 
When 



p[ln(p)-l] + l 



5" (i/, 1) — > 0, and ()118j) is not defined there. We consider asymptotic solutions of (jHl) of the 
form 



F,{x)=Ff\y) 



1 + A 



(119) 



where j = N — k as before. Using ()119|1 in (jH)) gives the following equation for Kf'\y) 

Kf - Kfl, = 

with solution 



Hence, 



F^iy) 



1 + a; j\ 



X-^k^^\y). 



To find k^^\y) we shall match p2()j) with the transition layer solution ()1()2|1 . 
As z ^ 1, (fin^ becomes 



[1^ (y, 2) , z] ~ y/eK{z) exp 



/27r 



exp ( --r 



2 'dr 



(120) 



(121) 



with 



From ()12Up we get as j 



V{y,l) 



y-Yo{l) 1 
1^2(1) x/^' 



1-2 



00 



Fi^\y) ~ v^/«(-2)exp 



-$(2) 



k^'\y). 



(122) 
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Matching ^TT^ and (fT^ gives 



V{y,l) 



27r 



and we conclude that 



^ ^iV\^' 1 1 



V{y,l) 



, , , , , exp ( — 1 dr 



(123) 



where 



[1 - 7) (p - 4A + 1) , 2A-7 + (7-l)A^ j^^^^ 



(A+ir 



(A + i; 



6.4 The boundary layer at z = 1, ^o(l) < y < oo 

To find the solution valid in the region 1 — z = 0{e), Yq{1) < y, we immediately conclude 
that Fk(x) must be of the form 



A \^/iV 



1 + A 



X-'-Fj'\y) 



(124) 



where Fj'^\y) is given by (|118p . This solution matches to and (jl23p as y Yo{l) and 
as j oo, respectively. 



7 The boundary x = 



For X = and A; < [cj , the values of Fk{0) can be computed from the ray expansion, since 
Ffc(O) - Ffc(oo) ~ £:K(0, z) exp [^^^(0, 2)] is well defined. For x = and k > [c] + 1, we have 
-^fc(O) = by (0}. We now examine how this boundary condition is satisfied by considering 
the scale x = 0(1) {y = 0(e)) and constructing a boundary layer correction to the ray 
expansion. Note that this part of the boundary is in the region i?*-^. 



7.1 The boundary layer at a: = 0, ^ < z < 1 

We shall find the solution satisfying the boundary condition ((7j). This boundary condition 
must be applied on the original x-scale. From ()93p we find that ^{y, z) = \l/*^'')(t/, z) + o{y), 
as 1/ — 0, where 



: (2; — 7) In 
7 ln(7) + 



ye 



i^-l)' 
^In 



+ (z - 1) ln(l - 2) - In (A + 1) + z\n{\) 



(125) 



Z — 'J 



ly 



+ 



z — 7 



[A(l - 27) + (A - 1)^] 
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Hence, we shall consider asymptotic solutions of the form 



Using ()12(ij) in ^ and taking into account that y = ex we get, to leading order. 



+ ^ — + 



2 + ;z - 37 
2{l -z){z-^) 



dz z — 'J dx 
The most general solution to (jl27|) is 



K^'^ = 0. 



x 



2; — 7 



Hence, 



i/7 



X, 2j ~ e ' exp 



-^'^'^\ex,z) 



1(1_^)IA;(7) (h). 



(126) 

(127) 

(128) 
(129) 



To find A;'-''^ (H) and z/y we will match ()129p with the corner layer solution ()71|). 

Recalling that / — a = and using the asymptotic formula for the Bessel functions PP , 



we get as e — s> and 'd fixed 



1 



2% - 7 



exp 



1 /e2; 



V27n7 V2z/ 



In 



2^ {z - 7) 



7? 



(130) 



Using (jl3(Jj) and writing (f7T|) in terms of x = and 2; = 7 + (Z — q;)£:, we have 



P 



27r Y 07(1 — 7) (2; — 7) 



X exp 



2; — 7 



A(l -7)e£ 
(z-7) 



In 



X exp < — In 

V 



exp < - 



1. 



7ln( ^ ) -(l-7)ln(l-7)-ln(A + l) 



2TTiJ ^ \^ 

Br 



— / -r + 1 - a exp - [x^? - (2; - 7) In (^9)] 



(131) 



7£ 



x^-7). 



a In - + 



2A(l-7) 
?9 



To evaluate (jl31|) asymptotically as e we shall use the saddle point method. We 
find that the integrand has a saddle point at = i, so that 



271/ y 07(1 — 7) 



exp <^ - In (A) + (^ - 1) ln(l - 7) - ln(A + 1) - 7 ln(7)] 



X exp 



2r — 7 



In 



2; — 7 



+ (j)E In 



75 



0(2; -7) 



(132) 



+ a In (0S) + 2A (1 - 7) H ^ -ET (0S + 1 - a) 

X 
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Taking the limit in ()129p as x — > 0, 2; — » 7 with S fixed, we obtain 



2; — 7 



(133) 



X exp 



z — 7 



In 



z — 7 



+ - [(z - 1) ln(l - 7) - In (A + 1) + z ln(A) - 7 ln(7)] 



Matching (fHI^ with (fH!^ we have 
and z/7 = |. Therefore, for 7 < 2; < 1, 



(1-7^ 



■ exp [{a - 0) S In + 05] 



e \2 
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X s 



P 



1 



07 (1 — z) 2; — 7 \z — 7/ \2; — 7 



r 



+ 1 - a 



1 



X exp - (2; — 7) In 



ee 



X exp 



(px 
2; — 7 



In 



7£ 



[{z - 1) ln(l - ;z) - In (A + 1) + ;2 ln(A) - 7 ln(7)] 

(134) 



0(2; -7) 



X 



+ 2A(1 -7) + (A - l)x 

2; — 7 



Note that from (IT!^ we have Fk{x) = O (x'^-L^J) , as x ^ 0, > [cj + 1. 

7.2 Matching the boundary layer at ^ = 0, ^ < z < 1 and the 
solution 

Writing x = j, and using Stirhng's formula we have, as e — > 

V2^ 



r 



Hence, 



(2; - 7) e 

X exp 

becomes, for x = - - 



+ 1 - a 



e (^ - 7) 



In 



(2 - 7) e 



e (2 — 7) £ 



(2; -7) 



P 



7(1-2) (z--f) 



X exp - (z — 7) In 



X exp 



£ (2; - 7) 



In 



2 
ye 

(z-7)2_ 



+ - [(2 - 1) ln(l - 2) - In (A + 1) + 2 ln(A) - 7 ln(7)] 



(135) 



+ 



A - 1 - 



z — 7 
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From 



and we get, as ?/ — > 



z — 7 
s ~ h 



In 



Using (IT!??)|) in ^ and (jHH) we find that 

1 



+ [A(2 - 37) + 7 + (A - l)z\ 



(136) 



K{s,t) 



P 



277 Y 7(1 - Z) 

in perfect agreement with p35|) . 



7.3 The corner layer at (0, 1) 

For {y, z) close to (0, 1), we use the variables x = - and j = N — k, and look for asymptotic 
solutions of the form 



:^(8)( 



Kf\x) 



with 



^!^^\x]e) = (1 -7) In 



xee 



.(1-7)^ 

Using ()137j) in gives, to leading order. 



In (A+ 1) + ln(A) -7ln(7). 



xX{3 + l)Kf^, = {l-^fKf 



whose solution is 



Kf\x) 



(i-7r 

Ax 



-M'^ ( 



X 



Hence, 



Fj^\x) ~ exp 



As j — > 00 ()138|) gives, by Stirling's formula. 



exp 



Ax 



U8-2j 



jyA:(«) (x) 



(i-7r 

Axj 



X 



(137) 



(138) 



(139) 



We determine k^^^ (x) by matching p39|) to the boundary layer expansion in ()134|) . 
Writing z = 1 ~ je and letting 2; — 1, we obtain from ()134j] 



-2i 



(1-7)^ 
Axj 



(140) 



X exp <! -^(^)(x; e) + In 
e 1 — 7 



£7 



0(1-7)J 



+ x(3A- 1) 
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By comparing ()139p and p40|) we find that 



x] 



r I \- 1 — a 



271 VI -7 
and z/g = 1. Therefore, 




pi \ <px 

exp < in 



~ — exp j -^(^)(a;; e) + In 
27r I £ 1 — 7 



1-7 



75 



1-7 



£7 



0(1-7), 



+ a;(3A- 1) 



X r 



(f>X 

1^7 



+ 1 - a 



(l-7)0j 

0X \ ° / p 1 

1 — 7/ Y07l — 7 



+ (3A- l)x 



(i-7r 

Ax 



-2i 



We can also show that the above, when expanded for x — ^ cxd, matches to the boundary layer 
expansion in ()118|) . valid for j = 0(1) and < y < Yq{1). 



8 The marginal distribution 

We will now find the equilibrium probability that the buffer content exceeds x 



N 



Pr [X{oo) >x]= M{x) = 1 - ^Ffc(x) 



(141) 



fc=0 



for various ranges of x. We will compare our results with those obtained previously by 
Morrison |31j . 

8.1 Approximation for x = 0(e) = 0{1/N) 

In this region we shall use the spectral representation of the corner layer solution which 
applies for x = ex = 0{e). Using the generating function 



J = -OC, 



xjz-' = exp 



X / _ 1 

2 r ~ i 



we obtain 



exp 



-(i+i) 



/=— oo 



4(j+i-«) 



Uo) . 
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Therefore, 

M{x) = M^^\x) ~ ^flj-exp 



X exp 



e. —i^il) exp 



X exp 



-$(7) 



j + 1 - a 
4^ j! (j + l-a) 



j + l-a' 



X exp 



^(j + l-a) + T(^,0 



Mo) 



From we have 

exp[T(^9,)] = [- (j + l-a)f+^xp 

Hence, 

M«(x)~v^,/?'^(7)exp 



2A(l-7) 



(j + 1 - a) 



1<|.(7) 



X > ^exp 



j>0 







j + l-a' 



X exp 



(142) 



The formula ()142|) agrees with Morrison's result (4.14) in [3T], taking into account the fol- 
lowing notational equivalences 



Morrison 


Dominici-Knessl 




1 — a 


r 


p — q 
7(1-7) ^0 


/t(7) 


_$(^) 


/(7) 


-ln(wo)-2£ 



8.2 Approximation for x = 0{£ ^) = 0{N) 

We shall now use the asymptotic solution in the region i?, as given by (jHS))- We have 



M(x) = M^'\y) = -J2G [y, ~ -J exp 

fc=0 n 



K(?/, z)(iz 



and using the saddle point method we get 



M(2)(y) 



V-^'..(y,7) 
37 



exp 



:^(2/,7) 



K(2/,7)- 



We recall that for a fixed y, z) is maximal at 2; = 7. From ()38|) and ()39|) we get 



where ?/s is understood to be evaluated at s = S{y,j) < 6*0 < 0. Thus 

From (1211) and dHHl) we get 

2/ (s, ^7) = 
and from (|93p and (|39p we have 



(143) 



+ p(A + 1)] + 2p 
A2 



T,) = sy {s, T,) - ln(A + 1) + ^ [(27 - 1) s - (A + 1)] + ^ In 



As 



p + S7(l -7) 



The results above agree with Morrison's (5.15) in ^SJm if we reconcile notation as below 



Morrison 


Dominici-Knessl 


r 


— s 


Z{r) 




In [Y{a)] 


A{-s)T,{-s) 


U{t) 


sy (s,T^) - ^(s,T^). 



9 Summary and discussion 

In most of the strip D = {{y, z) : y > 0, < 2; < 1} , the asymptotic expansion of Fk{x) 
G{y,z) is given by 



G{y,z) ~ eexp 



-'^{y,z) 



K{y, z) in R 



c 



or 



G(oo, z) — G{y, z) ~ —e exp 



-^{y,z) 



K(y, z) in R. 



(144) 
(145) 



If we consider the continuous part of the density, given by 

dG 

fk{x) = Fl{x) = e—{y, z), x > 0, 

oy 



38 




the transition between R and R'" disappears, and we have 



fk{x) ~ e'^y{y,z)e^^ 



e exp 



sKis,t), 



everywhere in the interior of D. Note that K{y,z) becomes infinite along y = Yq{z) (i.e., 
s = 0), but the product "^yiy, z)K{y, z) remains finite. 

The asymptotic expansion of the boundary probabihties -Ffc(O), k < [cj can be obtained 
by setting ?/ = in ()145|) . This expression can be used to estimate the difference 



Ffc(oo) - Ffc(O) = Pr X(oo) > 0, Z{oo) = k 



z < 7 



which is exponentially small for £ — * 0. Also, for a fixed z G [0, 7), fk{x) is maximal at x = 
(see Figure 4). 

However, for a fixed z G (7, 1), fk{x) is peeked along the curve y = Yq{z) (see Figure 5). 

This means that given k = zN > c active sources, the most likely value of the buffer will 
he X = NYq{z). If zN < c, the buffer will most likely be empty. For a fixed a; > 0, fk{x) 
achieves its maximum at 2; = 7 (see Figure 6). 



39 




40 
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Below we summarize the various boundary, corner and transition layer corrections to the 
results in fl92|) and ()91|) . where the paragraph number refers to the corresponding region (see 
Figure 7). 



1. k = I + c — a, X = ex, X = 0(1) 



l-a 



e 



-$(7) 



x^l e^'^ir it + l-a\ J,_^^^ f ^ ) exp [T(^)] dd. 



2TiiJ ^ \^ 

Br 



where J.(-) denotes the Bessel function, Br is a vertical contour in the complex plane 
with Re(^9) > and 



a = c— c 



7 + A — 7A, p = 7 — A + A7 



$(2) = -z \n{z) -{1- z) ln(l - z) + z ln(A) - ln(A + 1) 



Wo = -(1-7), = ]— } 

7 V27r - z) 



T(^9) = ( ^ - a 1 In 



2A(l-7) 



2-0 



In (mo) . 



2. y-ro(^) = 0(v/i), 7<^<1 

Ffc(a;) = ^) ~ v^/t(z) exp 



V 



U{z) 



with 



y - Yoiz) 



eJYoiz) 



I \ ^ " 7 P f ZX + Z - X\ 



A + 1 (A + 1) 



P 



z — 7 



(A + l)(Az + z-A)^ 



(A + lf V P 

+vT^(^-7) + (A-l)(.-7)^' 



L(A + lf (A + 1) 
( = 2A-7 + (7-l)A' 
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Yo(l) 




m 

y 



Figure 7: A sketch of the different asymptotic regions. 
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3. k^O{l) 

Fk{x) - Ffe(oo) = Flf\y) - Ffe(oo) - s^-^^exp 



-*(y,o) 



[A - 7^(2/, 0)]'^ - 



A- 7^(1/, 0) 



v^5(i/, 0) V [7(1 - 'y)S{y, 0) + p] Jo{y) ' 

4. A; = 7V-j, j = 0<y<yo(l) 



l + (l-7)>g(y,l) l^' 1 
A J j! 



X 



[l + (l-7)5(y,l)] 



^5(2/, 1) V [P + 7 (1 - 7) -5(2/, 1)] My) ■ 



A 



TV 



iVA-' 1 1 



V 



6. k = N-j, J = 0(1), y>Yo{l) 



exp ( -^T^ ) dr. 



A 



1 + A 



A-^-~Ff (1/), 



where F^^^ (y) is as in item 4 
7. x = 0{l), -f<z<l 

F,{x)^F^'\x,z)^[^f 



P 



(f)X 



(f)X 



+ 1 - a 



X exp - (2; — 7) In 



07 (1 — z — 7 \;s — 7/ \^ — 7 
""""^ V ^ [(^ - 1) ln(l - ^) - In (A + 1) + ^ ln(A) - 7 ln(7)] | 



X exp 



(j)X 

z — J 



In 



(z-7)2 

7£ 



- 7) 



+ 2A (1 - 7) + (A - l)x 



2; — 7 



8. A; = 7V-j, j = 0(l), x = 0(l) 



-l-2j 



Fit(x) = Ff (x) ~ ^exp <^ + -^In 
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X 



(i-7r 

Ax 



(f)X 



+ l-a 



(px 



(1-7)0. 

^ 1 

07 1 — 7 



+ (3A-l)a:^i 
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with 



*(^)(x;£) = (1-7) In 



exe 



(i-7r 



-7ln(7) +ln (^j^ 
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